Relaxation processes in a disordered Luttinger liquid 
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The Luttinger liquid model, which describes interacting electrons in a single-channel quantum 
wire, is completely integrable in the absence of disorder and as such does not exhibit any relaxation 
to equilibrium. We consider relaxation processes induced by inelastic electron-electron interactions 
in a disordered Luttinger liquid, focusing on the equilibration rate and its essential differences from 
the electron-electron scattering rate as well as the rate of phase relaxation. In the first part of the 
paper, we review the basic concepts in the disordered Luttinger liquid at equilibrium. These include 
the elastic renormalization, dephasing, and interference-induced localization. In the second part, we 
formulate a conceptually important framework for systematically studying the nonequilibrium prop- 
erties of the strongly correlated (non-Fermi) Luttinger liquid. We derive a coupled set of kinetic 
equations for the fermionic and bosonic distribution functions that describe the evolution of the 
nonequilibrium Luttinger liquid. Remarkably, the energy equilibration rate in the conducting disor- 
dered quantum wire (at sufficiently high temperature, when the localization effects are suppressed 
by dephasing) is shown to be of the order of the rate of elastic scattering off disorder, independent 
of the interaction constant and temperature. 

PACS numbers: 71.10.Pm, 73.21.-b, 73.63.-b, 73.20.Jc 



I. INTRODUCTION 

Strongly correlated electron systems in one dimension 
(Id) have become an area of immense interest from the 
perspective of both fundamental and technological as- 
pects of nanophysics. Intense experimental effort has 
focused on such realizations of quantum wires with a 
few or single conducting channels as cleaved-edge V- 
groove Q , and crystallized-in-a-matrix 0] semiconductor 
quantum wires, coupled quantum Hall edges running in 
opposite directions [4|, |aL single-wall carbon nanotubes 
p, polymer nanofibers [3], and metallic nanowires 
Central to much of the fascinating physics of the Id sys- 
tems is that electron-electron (e-e) interactions in Id ge- 
ometry can have dramatic effects leading to the emer- 
gence of a Luttinger liquid (LL) [l(| • The latter con- 
stitutes a canonical example of a non-Fermi liquid, in 
which quasiparticle fermionic excitations are inappropri- 
ate to describe low-energy physics. At the foundation of 
the conventional LL theory is a description in terms of 
bosonic elementary excitations (plasmons, spinons) [Io| . 
Following this approach, the ground-state properties of 
a clean LL are well understood for arbitrary strength of 
interaction. Much has also been learned about the LL in 
the presence of a single compact scatterer 1(| 11 1. How- 
ever, as far as a disordered LL is concerned, a number of 
important questions, even at the most fundamental level, 
remained largely unanswered until very recently (for a re- 
view see Ref. (l3|). 

In the presence of disorder, quantum interference of 
scattered electron waves leads to the effects of Anderson 



In a Id electron gas, all electron states are localized even 
for an arbitrarily weak random potential and the local- 
ization length is the mean free path. In the case of non- 
interacting electrons, the quantum localization in Id has 
been studied in great detail (see, e.g., Ref. [3). A princi- 
pal complication that arises in the disordered LL is that 
the quantum interference phenomena yielding the An- 
derson localization are conventionally treated in terms of 
fcrmions, by employing the concepts of interference and 
dephasing of fermionic excitations. The question of to 
what extent the notion of phase relaxation in the local- 
ization problem is applicable to the (non-Fermi) LL is 
therefore of crucial importance. Recently, this problem 
was addressed in Refs. 12, llal, where the interaction- 



localization 13[. Similarly to e-e interactions, the lower 
the dimensionality, the stronger the localization effects. 



induced dephasing rate that governs the localization term 
in the conductivity of the disordered LL was calculated. 

Another conceptually nontrivial aspect of the interplay 
between disorder and interaction concerns the nonequi- 
librium properties of the LL. In the homogeneous case, 
the LL model is completely integrable and as such does 
not exhibit any relaxation to equilibrium: an arbitrary 
excited state will never decay to the equilibrium state 
characterized by temperature. Of central importance 
is therefore the question of how the equilibration of 
fermionic and/or bosonic excitations in the LL occurs 
in the presence of disorder. 

This paper is primarily concerned with various relax- 
ation processes associated with inelastic interactions be- 
tween electrons in the disordered LL. Specifically, we fo- 
cus on the rates of e-e scattering, phase relaxation, and 
energy relaxation, with emphasis on the essential differ- 
ences between them. In Sec. II we begin with the formu- 
lation of the model. Section III highlights a few aspects of 
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temperature-dependent screening of disorder in Id. Sec- 
tion IV covers the problem of phase relaxation — this 
discussion largely follows the results of Refs. 12, HI and 
serves as the starting point for our approach to nonequi- 
librium physics of the LL. In Sec. V we consider energy 
relaxation and introduce a general framework [l6| for 
studying the behavior of the disordered LL out of equi- 
librium. 



II. MODEL 

Let us specify the model. By decomposing the elec- 
tron operator into right- and left-moving parts , ip{ x ) = 
ip+(x) + ip-(x), the Hamiltonian of a disordered LL is 
written as 



H 

Hkin 



-Vf ^2 dx ^i (i(J>dx + k F ) Vv i 

\j2 J dx ( n » v f + n ^ n M 

\u b (x) ^>++H, 



dx 
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(4) 



Here n M = t/'u^/j is the density of the right and left 
movers and their dispersion relation is linearized about 
two Fermi points at the wavevectors ikp with the ve- 
locity v F - Throughout the paper we consider spinless 
electrons (for spin-related effects see Ref. |17|). 

The e-e interaction, Eq. (|3|), is characterized by the 
Fourier components of the short-range (screened) interac- 
tion potential with zero momentum transfer Vf (forward 
scattering between right and left movers) and Vf (for- 
ward scattering of electrons from the same chiral branch 
on each other). Unless the right and left movers are 
spatially separated (as in coupled quantum Hall edges), 
Vf = Vf. The Luttinger model per se does not include 
backward scattering characterized by the Fourier com- 
ponent Vb with momentum transfer ±2fcp. For spinless 
electrons, however, Vb can be trivially incorporated by 
shifting Vf — > Vf — Vb, since two types of scattering — due 
to Vf and Vb — are then related to each other as direct 
and exchange processes. The local interaction between 
identical fermions Vf yields no scattering, but, due to a 
quantum anomaly in the LL model, generates a shift of 
the Fermi velocity vp — » v* F = v F + Vf/2-K. It is cus- 
tomary to parametrize the strength of e-e interaction by 
means of the Luttinger constant K: 

A- = (-^| 1/2 . n (5) 
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The velocity of elementary excitations (plasmons) in a 
clean LL is given by 



which transforms for Vf = Vf into u = vp/K, 

The low-energy theory described by the Hamiltonian 
(fTJ) is only then well-defined when supplemented by an 
ultraviolet energy cutoff A. The latter depends on mi- 
croscopic details of the problem and obeys 



A = u/wX 



(7) 



where the length scale A is set by the lattice constant, 
the Fermi wavelength, or the spatial range of interaction 
in the original microscopic theory, whichever gives the 
smallest A. Thus the complete set of parameters defin- 
ing the LL model in the absence of disorder includes v* F , 
Vf, and A. It is worth noting that the input parame- 
ters of the low-energy theory include Fermi-liquid-type 
renormalizations coming from energy scales larger than 
A; in particular, the "bare" vf in Eq. |T]) in general is not 
an interaction-independent constant if the interaction is 
strong (I - K ~ 1). 

The term ffdis, Eq. ([4|, describes backscattering of 
electrons off a static random potential U(x). The lat- 
ter is taken to be of white-noise type with the correlators 
of backscattering amplitudes 



U b (x)U^{0) = U(x)U{0) = wS(x) 



(8) 



and Ub(x)Ub(0) — 0. The forward-scattering amplitudes 
are omitted in Eq. ((¥]) since they can be gauged out in 
the calculation of the conductivity. 



III. ELASTIC SCATTERING 

One of the characteristic features of a LL is a large 
renormalization of the strength of disorder ([5]) by e-e in- 
teraction. In particular, the conductivity without any 
localization [l3| or pinning [l8| effects included ("Drude 
conductivity") is od{uj,T) = e 2 VF/ir{—iuJ + M(uj,T)], 
where the disorder-induced scattering rate in the dc limit 



f 



r(T) 



ReM(0, T) 
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T 



2(1- K) 



(9) 



grows as a power law with decreasing temperature T 
for repulsive interaction (K < 1). The momentum re- 
laxation rate in the absence of interaction is given by 



2wv F with w from Eq. 



Calculating the 



^(l-a 2 ) 1/2 



(6) 



Drude conductivity at finite oj and sending uj — > after- 
wards allows to unambiguously determine the coefficient 
8^ = ^(1 + K)/Y(2K ) in the relaxation rate. Here 
and below the disorder is supposed to be weak in the 
sense that At ^> 1. 

The underlying physics of the renormalization ([9]) can 
be described in terms of the T-dependent screening of 
individual impurities; specifically, in terms of scattering 
by Friedel oscillations which slowly decay in real space 
and are cut off on the spatial scale of the thermal length. 
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At this level, the only peculiarity of the LL as compared 
to higher dimensionalities is that the renormalization of 
r is more singular and, most importantly from the cal- 
culational point of view, necessitates going beyond the 
Hartree-Fock approximation, even for weak interaction 
(see, e.g., Ref. [20(). 

In general, not only the strength of disorder but also 
the strength of interaction is subject to renormalization 
and depends on T, so that the function r(T) is not a 
simple power law. An important question, therefore, is 
under what condition the exponent in Eq. @ is given 
by the bare interaction coupling constant. One of the 
approaches to the problem was formulated in Ref. [2lj 
in terms of a bosonic renormalization group (RG). The 
RG approach does not allow to obtain the A-dependent 
prefactor olk in Eq. ([5]), but is particularly beneficial in 
predicting the T dependence of the Drude conductivity. 
For spinless electrons, the one-loop RG equations read 

dK/dC = f{K)V , (10) 
dV/dC = (3 - 2K)V , (11) 

where C = InL/A and V = 2wX/ttu 2 . For the Drudc 
conductivity (i.e., as long as the localization effects are 
not included, see Sec. IV), the spatial scale L is given 
by the thermal length u/T. The scattering rate l/r(T) 
is then proportional to TT>(T). The function /(A) — 
-A 2 /2 + (l + A 2 )(3-2A)/4 vanishes at A = 1, so that 
interaction is not generated by disorder (in the original 
equations of Ref. [2lJ , the coupling constant A contains 
an admixture of disorder and therefore the corresponding 
/(l) 7^ 0, see Ref. [H} for a discussion of this point); 
moreover, the interaction (hence 1 — A") does not change 
sign in the course of renormalization. 

The RG flow pop .lfTT jl is characterized by a separatrix 
which behaves as V = 8(A - 3/2) 2 /9 for A > 3/2 and 
terminates at A = 3/2. For the bare (taken at L = A) 
values of T> and K that lie below the separatrix (i.e., for 
the case of strong attractive interaction with A > 3/2), 
the disorder strength T> renormalizes to zero, otherwise T> 
grows with increasing L to a strong-coupling point with 
T> ~ 1. The renormalization of the coupling constant 
A by disorder is essential if the RG trajectory is close to 
the parabola V = 8(A-3/2) 2 /9. For example, if the RG 
flow passes through the point A = 3/2, the integration 
of Eqs. (HOJ and (HJ gives for V < 1: 

/3-D 1 / 2 l\ 
V-V =V tan 2 I In— , (12) 

where Vq and Lq are the values of V and L at A = 3/2 
and the sign of ln(L/L ) is positive for running A < 
3/2 and negative otherwise. One sees that T> grows with 
increasing L for A < 3/2 as 

V = 8/9\n 2 (l/L) (13) 



(for T>q <C T> <C 1). Here the renormalized mean free 
path I (the scale at which T> ~ 1) obeys \sx{1/Lq) = 
2 1 / 2 7r/32?y 2 . The logarithmic dependence of 2? on i 
is precisely due to the renormalization of A. 

On the other hand, if the bare A < 3/2, the RG tra- 
jectory follows Eq. (fT2|) with Lo = A and Vq understood 
as the bare value of V only at V - D > (AT - 3/2) 2 . 
Integrating Eqs. (fTQ)) . (|lip in the opposite limit 

V - V < (A - 3/2) 2 , (14) 

one gets 

V = V {L/\f~ 2K , (15) 

which corresponds to Eq. ^j. Equation (Ti~4"| thus answers 
the question of when the renormalization of A may be ne- 
glected. Notice that for repulsive interaction (A < 1) the 
condition (fT4|) is satisfied for the whole range of T> <C 1 
[which is where the RG equations dTOj) , (jTTJ) are valid]. It 
follows that for the most relevant case of direct Coulomb 
interaction the renormalization of interaction on ballistic 
scales (C< 1) plays no role and the exponent in Eq. © 
is T-independent and given by the bare value of A (the 
one in a clean system). In other words, the renormaliza- 
tion of disorder for repulsive interaction reduces to the 
renormalization of an individual impurity. It is worth 
emphasizing that this does not mean that the disorder- 
induced correction to the bare value of 1 — A is small: in 
fact, the correction is of the order of 1 — A itself when 
T> ~ 1. The point is that the exponent of T>{L) and, cor- 
respondingly of the renormalized scattering rate 1/t(T) 
is not given by the running coupling constant A, but 
rather is accumulated on the whole RG trajectory. 

IV. PHASE RELAXATION 

The renormalization of r stops with decreasing T at 

Tr(T) ~ 1 , (16) 

since the long-range Friedel oscillations created by dis- 
order are cut off even at zero T on the spatial scale 
of the disorder-induced mean free path. This condi- 
tion gives the zero-T mean free path I cx Tq^ 3 2K ^ [no- 
tice that Eq. (fTB)) is also expressible as T>(L) ~ 1 with 
L = u/T] and, correspondingly, the zero-T localization 
length £ ~ I. It is important to stress, however, that 
the above condition does not correctly predict the on- 
set of localization with decreasing T — in contrast to the 
argument, frequently stated in the literature (see, e.g., 
Ref. and references therein) and based on the RG 
equations (fT3)) . (|14p . which treat scalings with the length 
scales L and u/T as interchangeable. While substituting 
u/T for L is justified for the "elastic renormalization" 
[Eq. ©], the one-loop equations (fTB")) . (TT4")) miss, by con- 
struction, the interference effects (coherent scattering on 



4 



several impurities) that lead to localization. The status 
of the RG (2l| is thus that of the Drude formula for in- 
teracting electrons. The T dependence of the conductiv- 
ity tr(T), however, comes not only from the T-dependent 
screening of disorder [Eq. ©] , but also from the localiza- 
tion term in a(T) whose amplitude is governed by phase 
relaxation due to inelastic e-e scattering. The tempera- 
ture below which the localization effects become strong 
is, in contrast to Eq. (I16p . determined by the condition 



order in the interaction is given by 



-(r)/7*(r) ~ 1 



(17) 



where is the weak-localization dephasing time. Notice 
that for weak interaction (1 — K ~ a <C 1), Eq. (fT7|) 
is satisfied at much higher T than Eq. (fl6|) . Below we 
introduce the notion of dephasing of localization effects 
in the disordered LL and analyze the phase relaxation in 
the limit of weak interaction. 

The very applicability of the notion of dephasing, as 
we know it from the studies of higher-dimensional Fermi- 
liquid systems, to the LL is not altogether apparent. A 
subtle question concerns the nature of elementary exci- 
tations in the LL, especially in the presence of disorder. 
The clean LL is a completely integrable model which is 
represented in terms of noninteracting (hence nondecay- 
ing) bosons; however, the phase relaxation in electron 
systems is conventionally described in terms of inter- 
acting fcrmions. Physically, the difficulty is related to 
the fact that the bosonized approach describes propa- 
gation of density fluctuations, whereas the natural lan- 
guage for quantum interference phenomena is that of 
quantum amplitudes. To study the interference effects 
and their dephasing, one has therefore to either proceed 
with the standard bosonization, poorly suited to describe 
the quantum interference in the inhomogeneous case, or 
try to define the observables in such a way that they 
can be expressible in terms of decaying fermionic excita- 
tions. In what follows in this section, we take the lat- 
ter path and give a succinct analysis of the phase relax- 
ation in the disordered LL, based on the results obtained 
within the "functional bosonization" formalism [12| and 
the quasiclassical formalism 115], both of which combine 
the fermionic and bosonic approaches to the problem. 

Let us first point out one of the subtleties of the LL 
model, which is crucial to our discussion of the phase 
and energy relaxation. The Golden rule expression for 
the e-e collision rate at equilibrium, as follows from the 
Boltzmann kinetic equation, reads 



dui / de' K{ui) 



Tee(e) 



(fe-u)fe'fe'+u)+fe-uf e 'fe'+u) , (18) 



where / e is the Fermi distribution function and ff? — 
1 — f e . Consider the clean case. Then the scattering 
kernel K(u) = /Cf + (w) + fCf_{u) + K, f {uj) to second 



7T 3 /9 



dq 
2^ 



ReD + (io,q)[ 



Vf f dq 

= ~r hr ReD +^'^ ReD -^' 

ir°p J Zir 



(19) 



cj,q), (20) 



and K, F = — /C+ + - The Hartree terms /C+ + and )C+_ arc 
related to scattering of two electrons from the same or 
different chiral spectral branches, respectively, K, F is the 
Fock (exchange) counterpart of /C+ + , the thermodynamic 
density of states p — l/irup, and D± — iirp/(oj^fqvp+iO) 
are the two-particle propagators in the clean limit. 

Substituting Eqs. Ifl9 |) ,(p0 || in Eq. JTSJ), we obtain the 
lowest-order result for the e-e scattering rate at the Fermi 
level (e = 0) in terms of the corresponding contribu- 
tions to the retarded electronic self-energy S + defined 
by G+(e,p) = [e — vpp — S + (e,p)]~ 1 , where G+ is the 
retarded Green's function for right-movers. Specifically, 
t- 1 = -2Im£+(0, 0) with E+(0, 0) = Sf + + Ef_ + E F , 
where 



T n 2 f i ( i w i w \ 

ImL , _i_ = a vf dujui [ coth tanh — 

+± 2 J V 2T 27V 

/ dqS(uj — vpq)5(u! =p vpq) , (21) 



S F = — and we put Vf = Vf . One sees that the 
contribution of is diverging. For spinless electrons, 
however, the divergency is canceled by the exchange in- 
teraction. Indeed, as we have discussed in Sec. II, the Vf 
interaction drops out of the problem in this case, induc- 
ing only a shift of the velocity vp — > v F . It is worth not- 
ing that the "Hartree- Fock cancellation" is only exact for 
the point-like interaction (when Vf is independent of the 
transferred momentum), otherwise Vf yields a nonzero 
contribution [22] to t^ 1 . The latter is small and can be 
neglected in the low-T limit for Vf = Vf but is the only 
one present for Vf = 0, which is the case, e.g., for an 
isolated quantum Hall edge. The remaining term 
gives 



T-} = -2ImE?_ = ira 2 T 



(22) 



This may look very similar to the familiar T 2 or T 2 In T 
dependence of the e-e scattering rate in clean three- or 
two-dimensional electron systems, respectively. However, 
the nontrivial point — which demonstrates the peculiarity 
of the LL model — is that r" 1 in Eq. (|2"2")1 is determined 
by 



u,q = o , 



(23) 



i.e., by scattering processes with infinitesimally small en- 
ergy transfers, in contrast to higher dimensions where the 
characteristic transfer is of order T. On the other hand, 
it is worth emphasizing that Tr ee 1 for a <C 1, which 
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in Fermi-liquid theory is commonly referred to as one of 
the conditions for the existence of a Fermi liquid. In this 
respect, the weakly interacting LL, while being a canon- 
ical example of a non-Fermi liquid, reveals the typical 
Fermi-liquid property. The LL physics (e.g., the power- 
law singularity of the tunneling density of states at the 
Fermi level) is in fact encoded in the singular real part of 
the self-energy £+(e,p) (for more details see Ref. fl2l|). 

It is the property (|23j) that actually makes the Id case 
special as far as the dephasing problem is concerned. In- 
deed, in the spirit of Ref. [23j, soft inelastic scattering 
with qvp,uj <C t^ 1 is not expected to contribute to the 
dephasing of the localization effects. In higher dimen- 
sions, in the ballistic regime TV 1, this infrared cutoff 

is given 



is of no importance and the dephasing rate t , 1 ; 



[24( by Tgg 1 . However, in view of Eq. ([231), r^ 1 in Id 
cannot possibly reduce to r~ . 

The dephasing rate can be accurately defined by 
calculating the weak-localization correction to the con- 
ductivity of the disordered LL as a function ofT [II El- 
The leading localization correction Act in the ballistic 
limit t^/t -C 1 is related to the interference of elec- 
trons scattered by three impurities. One of the diagrams 
contributing to Act (for the complete set of diagrams at 
the leading order see Ref. [H, 1H| ) is given by a "three- 
impurity Cooperon" (Fig. [1]), which describes the prop- 
agation of two electron waves along the path connecting 
three impurities ( "minimal loop" ) in time-reversed direc- 
tions. In the absence of dephasing, quantum interference 
processes involving a larger number of impurities sum 
up to exactly cancel (similarly to the noninteracting case 
[1J|) the Drude conductivity ctd = g 2 vpt/tt, where r is 
given by Eq. ©. For r^/r C 1, they only yield sublead- 
ing corrections through a systematic expansion in powers 
of t^/t. 

Within the functional-bosonization description of the 
LL [25|, extended in Ref. [HJ to treat disordered 
systems, the interaction can be exactly accounted for 
by performing a local gauge transformation ip^ix^r) — » 
ipnix, t) exp[i^(a;, r) ], where the bosonic field 0^(x,t) 
is related to the Hubbard-Stratonovich decoupling field 
(p(x,r) by 



(24) 



Here r is the Matsubara time. The correlator 
{<p{x, r)ip(0, 0)) = V{x, t) gives the dynamically screened 
interaction, for which the random-phase approximation 
(RPA) in the LL model is exact [27[. In the presence of 
impurities, the interaction can thus be completely gauged 
out to the backscattering impurity vertices — Eq. ([24)) is 
then exact for any given realization of the impurity po- 
tential. In Fig. [TJ the fluctuating disorder-induced gauge 
factors are denoted by the wavy lines attached to the 
backscattering vertices: each impurity vertex contributes 
the factor exp[±(# + — #_)] and the averaging over fluc- 



tuations of if pairs all the fields 0± with each other. The 
interaction thus induces the factor 



exp(-Sc) = (exp[i(6 f -e b )}) 



(25) 



to the Cooperon loop, where Of and Ob are the phases 
accumulated by an electron propagating along the "for- 
ward" and "backward" paths and the averaging is per- 
formed over the fluctuations of the field (p. Notice 
that the averaging couples with each other not only the 
phases 0± attached to the impurities shown in Fig. Q] 
but also those attached to impurities which yield damp- 
ing of the d yna mically screened interaction. As shown 
liL IH, the boson damping is crucially impor- 



Refs. 



tant for the dephasing (see below) and a parametrically 
accurate approximation is to include impurity-induced 
backscattering in the effective interaction at the level 
of the disorder-dressed RPA ("dirty RPA"). The total 
Cooperon action 



Sc — S + S IC 



(26) 



accounts then for both the dephasing (S) and the elastic 
renormalization of impurities (S ranorm ) and we refer the 
reader for technical details of the formalism to Ref. [l| • 
The leading localization correction to the conductivity 
can be represented in the form 13, HI 



Act 



-2ct d / dt c / dt a P 2 (t c ,t a )exp\ 
10 Jo 



S(t c ,t a )] 



(27) 

where P 2 {t c ,t a ) = (l/8r 2 ) exp(-< c /2r)9(t c - 2t a ) is the 
probability density of return to point x — after two 
reflections at points x = ut a and x = —u(t c /2 — t a ). 
Here ut c gives the total length of the Cooperon loop and 
ut a , being the distance between two rightmost impuri- 
ties, parametrizes the geometry of the loop. The classical 
trajectory for the Cooperon is characterized by a single 
velocity 12j and this is u (the difference between u and 
Vf can be ignored for a<Cl, but uniformly on the whole 
trajectory). The phase relaxation is encoded in the de- 
phasing action S in Eq. (|27[) , which is a growing function 
of the size of the Cooperon loop and cuts off the local- 
ization correction at t c ~ t^. The dephasing rate t^ 1 is 
thus defined by the characteristic scale of t c on which the 
dephasing action S ~ 1. 

In the limit t c t the action reads [13, ll_5j : 



S(t c ,t a ) = 2ira 2 Tt a (t c -2t a )/r 



(28) 



Substitution of Eq. ((28]) into Eq. ([27]) gives 
1 /t,a2, t_ 1 

where 

t-7 1 = a(nT/r)^ . (30) 



Ao = -\av( T ^) 2 \n-^^-Ho?T) , (29) 

4 \ T J Ta, Ct z ± 
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FIG. 1: Three-impurity Cooperon diagram with interac- 
tion effects encoded in the fluctuating factors exp(±0 (J ) (de- 
noted by the wavy lines) attached to the backscattering ver- 
tices (marked by the crosses). The dashed lines connect the 
backscattering vertices belonging to the same impurity (e.g., 
1 and 1 refer to two backscatterings off impurity 1 at two 
different times). 



One sees that the phase relaxation in the disordered LL 
occurs on time scales much longer than the lifetime r ee : 



= (T ee T) 1/2 > T, 



(31) 



Note that r^ 1 vanishes in the clean limit, in contrast 
to the total e-e scattering rate — in agreement with the 
observation ([23]) and the basic fact [23j that scattering 
with energy transfers smaller than t^ 1 is not effective in 
dephasing the localization effects. 

The vanishing of the dephasing action at t _1 — > 
can be made more transparent from the technical point 
of view by representing S as a difference between the 
self-energy (5g + Sbb) and vertex (Sib + Shi) contribu- 
tions. Here the terms Sij , associated with an inelastic in- 
teraction between electrons propagating along the paths 
Xi(t) and Xj(t), where i,j = (f) and (b) stand for the 
"forward" and "backward" time-reversed paths in the 
Cooperon, are given by 



duj f dq [*' , f'« 1 „ , , 
2^7 dh J dt.-lraV^q) 

x exp{iq [xi(t{) — Xjfa)] — uo(t\ — 1%)} . (32) 



S^ — T 



eqnarray 

Equation (|32p is similar to that in higher dimensional- 
ities ("AAK action" 23J) with one subtle and important 
distinction. Because of the Hartree-Fock cancellation of 
the bare interaction Vf between electrons from the same 
chiral branch [recall the discussion after Eq. (|2Tj) ]. the 
dynamically screened retarded interaction V^(u>,q) ac- 
quires the indices [A, v denoting the direction of motion 



of the interacting electrons: fi = sgnij, v = sgnij. If 
one would keep both Vf and Vf processes in V(u>,q), the 
dephasing action in Id could not be written in the form 
of Eq. (|32[) — in contrast to higher dimensionalities, where 



Sij is given by Eq. {32} with the "full" V(u,q). 

Neglecting the disorder-induced damping of the dy- 
namically screened interaction yields 



Se = Sib = t c /2r e , 



(33) 



and the exact cancellation of the self-energy and vertex 
parts in the total dephasing action, hence the vanishing 
of t^ 1 (|30|) in the clean limit. It is thus only because of 
the small difference between Ss and Sib produced by the 
dressing of V^ v (w, q) by impurities ( "dirty RPA" ]2 15 1) 
that the dephasing action (|2"8|) is not zero. The charac- 



teristic energy transfer w in the processes that lead to the 
dephasing (i.e., contribute to the difference Ss — Sib) is 
much larger than r _1 [more accurately, to is spread over 
the range between tT 1 and r _1 , because of the logarith- 
mic factor in Eq. (|29[) ]. which justifies the expansion of 
S in powers of r _1 , while the condition Tr^ 3> 1 justifies 
the quasiclassical treatment of the electromagnetic fluc- 
tuations in Eq. (32|). Substituting Eq. P| into Eq. (fT7|) 
gives the temperature scale T\ ~ l/a 2 r below which the 
localization effects become strong (for the behavior of the 
conductivity at T < T x see Ref. [23]). Note that T\t » 1 
for weak interaction. 



V. ENERGY RELAXATION 

We now turn to the noncquilibrium properties of the 
disordered LL [16j . Here we are primarily interested in 
the equilibration rate at which an excited state relaxes 
to equilibrium (other aspects of the nonequilibrium relax- 
ation will be discussed elsewhere fl6i] ) . As mentioned in 
Sec. I, the integrability of the clean LL model precludes 
energy relaxation. The absence of inelastic scattering in 
the LL deserves additional comment. For scattering of 
electrons from different chiral branches on each other, 
the energy and momentum conservation laws for linear 
electronic dispersion lead to two equalities: to — vpq = 
and u> + vpq = 0. These combine to give u>, q = and 
thus no energy exchange [cf. Eq. For scattering of 

electrons of the same chirality //, the energy-momentum 
conservation laws give a single equation u> — ^,vpq = and 
at first glance the energy relaxation is allowed. More- 
over, the relaxation might seem to be very strong since 
the Golden-rule expression for the probability of scat- 
tering contains the delta function 5(ui — jMVp) squared. 
For the point-like interaction, the diverging Hartree and 
exchange terms cancel each other; however, for a finite- 
range interaction — despite the LL model being still com- 
pletely integrable — the cancellation is no longer exact. 
The energy relaxation, nonetheless, is absent in the LL 
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model for a generic shape of the interaction potential. 
The point is that beyond the Golden rule the diverging 
terms sum up to produce the dynamically screened in- 
teraction between electrons (exactly given by the RPA), 
which propagates with velocity u(q) ^ Dp. As a result, 
the probability of scattering contains a product of two 
delta functions 6(tu — fiVFq)6[uj — (iu(q)q\ with different 
velocities, which yields u>, q — for electrons from the 
same chiral branch as well. 

The energy relaxation is thus only present if one goes 
beyond the clean LL model. One possibility comes from 
three-particle scattering 29], which occurs for a nonzero 
range of interaction provided that the electronic disper- 
sion is nonlinear. The three-particle collision rate is small 
in the parameter T/ep *C 1, where ep is the Fermi en- 
ergy. Another possibility is to take into account impurity 
backscattering, which may lead to a much stronger mech- 
anism of energy relaxation. 

It is important that the nonequilibrium state of the 
LL in general cannot be described in terms of a single — 
either bosonic or fermionic — distribution function. The 
simplest example to illustrate this point is that of the 
clean LL in which the left and right movers, separately 
at equilibrium within themselves, are characterized by 
the Fermi distribution functions = fp(e — fi±) with 
different chemical potentials. Then the distribution func- 
tions (lu) of the left and right plasmon modes are con- 
structed as the convolutions of the fermion functions: 



N ± (u) = - [ deffil-f; 
u J 



(34) 



i.e., are given by the equilibrium Bose distribution, inde- 
pendent of fi±. This observation shows that the purely 
bosonic description of the clean LL at a finite bias volt- 
age is not complete. Such a "partial nonequilibrium" 
setup, in which the bosons are still at equilibrium, has 
been studied previously by employing the conventional 
bosonization (see, e.g., Ref. [30(). Furthermore, the 
nonequilibrium transport through a single impurity be- 
tween equilibrium leads shifted by the voltage fi- has 
been studied in Ref. [31|. However, the standard scheme 
of bosonization will break down if the nonequilibrium dis- 
tributions of the injected right- and left-moving electrons 
are not the Fermi distributions. 

The challenge is thus to formulate a theoretical frame- 
work to describe a genuinely nonequilibrium situation 
in which both the bosonic and fermionic excitations are 
out of equilibrium. It is worth stressing that in the in- 
homogeneous case the nonequilibrium distribution func- 
tions do not obey the simple local relation Q34p , since the 
distribution functions of plasmons and electrons evolve 
with different velocities (u and vp). Notice also that the 
necessity of introducing both the bosonic and fermionic 
distribution functions is not pec uliar to Id: for higher- 
dimensional systems see Ref. [32j . 

Our approach to nonequilibrium phenomena in the LL 



uses as a base the formalism of the "functional bosoniza- 
tion", developed in Ref. [l2| for the treatment of dis- 
ordered LL at equilibrium. A conceptually similar for- 
malism has been applied earlier for higher-dimensional 
disordered conductors in the study of the single-particle 
density of states out of equilibrium in Ref. [33j. The 
nonequilibrium tunneling density of states in the clean 
LL has been considered within the functional bosoniza- 
tion approach in Ref. [3] . Here we formulate the theory 
of the disordered LL out of equilibrium, which builds on 
the approaches of Refs. 35|, 3g] and Ref. 32 1, in terms of 
the effective nonequilibrium real-time action. To account 
for the e-e interaction, we introduce a dynamical field 
4>{x, f) which decouples the four-fermion term in the ac- 
tion by means of the conventional Hubbard-Stratonovich 
transformation. 

The central object of our theory is the quasiclassical 
Green's function g(x, t\, t^) for electrons in the LL, taken 
at coinciding spatial points (37l |: 



g(x,t 1 ,t 2 ) = lim 2iv F g(x,x',tx,t 2 ) 

x' — >x . 

- sign(x - x')5(ti - t 2 ) 



(35) 



This function, which is a 4 x 4 matrix in the Keldysh 
and channel (right/left) space, satisfies the Eilenberger 
equation [38j]: 



where 



iv F d x g + [ id t T z -H ,g] = , 



H = 4>T z + -(U b T+ +U* b T-) 



(36) 



(37) 



Here and throughout this section below, vf means 
the renormalized velocity v* F [see the discussion around 
Eq. (J5J)], so that the difference between u and vp is of 
order a 2 for small a. Equation (f3"6")l describes the mo- 
tion of an electron in the random potential character- 
ized by the backscattering amplitude Ub(x) [Eq. {8])] in 
the presence of the dynamic field <p(x,t) = diag (0 + , </>~), 
where <^(x, t) = (f>i(%, t) + a x <p2(x, f) and 0^ and <f>2 are 
the classical and quantum components of the Hubbard- 
Stratonovich field with chirality /z, respectively. The 
Pauli matrices t z and — t x ± ir y act in the chan- 
nel space. We also introduce the Pauli matrices cf x ,y,z 
that act in the Keldysh space. The Hamiltonian H (|3"7| 
is defined on the direct product of the time, Keldysh, 
and channel spaces. Accordingly, the commutator [ , ] in 
Eq. (l36|) is understood with respect to all three ("dis- 
cretized" time, Keldysh, and channel) indices. The op- 
erator dt acts as 9^ from the left and as (— 9* 2 ) from 
the right. For the case of linear electronic dispersion, as- 
sumed in the LL model, the Eilenberger equation |36|) 
is exact for any given realization of the backscattering 
amplitude Ub(x). 
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The next step is to average Eq. (|36|) over disorder 
this point we disregard the localization effects 



At 



12, m 



which limits the applicability of the subsequent deriva- 
tion to sufficiently high (effective) temperatures; specif- 
ically, for the length of the quantum wire larger than 
the mean free path to T ^> T\ ~ l/a 2 r (see the end of 
Sec. IV). Under this condition we can perform the dis- 
order averaging at the level of the self-consistent Born 
approximation, which gives 



id t ~4>+ — g~v,gv 
4r 







(38) 



for the disorder-averaged Green's function gi*. In what 
follows we only deal with the averaged propagators and 
therefore omit the bar for brevity. 

The Green's function g satisfies the normalization con- 
dition 



gog= iS(t! - t 2 ) , 



(39) 



where the dot denotes the convolution in all three (time, 
Keldysh, and channel) spaces. The constraint (|3"9"|) en- 
ables us to formulate the effective action that reproduces 
Eq. (|38|) as its saddle point in the form essentially com- 
bining the actions derived in Refs. [35| and (36l |: 



S[g,(f>,A] = -— Tr (id t - <i>)r z + v F A 
2vf L 
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- l Trg Q T- 1 d x T--!—Trg + g-. 
2 8vpT 

(40) 

The Green's function in Eq. (|40|) is represented as g — 
TgoT- 1 = diag(<7+,-<7"), where g = diag (<ftj~, -<7o~) 
corresponds to the saddle point of the action of the non- 
interacting problem and the unitary transformation T 
(diagonal in the channel space) parametrizes possible 
fluctuations around go [satisfying the constraint Q39|)] ; 
induced by fluctuations of the field <j)(x,t). To gen- 
erate the response functions in the Keldysh formalism 
[39I . |40| . we have also added the external-source term 
A(x,t) — a\(x,t) + a x a 2 (x,t). The trace operation in- 
cludes the summation over the Keldysh and space indices 
and the integration over time. The Keldysh partition 
function of the system can now be expressed as a func- 
tional integral over <j), 



Z[A] 



V<ft )2 (x,t) exp{iS[(f>,g, 



A] 



iTr 
2 



2ttvf 



(41) 



where g(x, t\, ti\ [(f>(x, t)]) is understood as minimizing 
the action (|40|) for a given <f>(x, t) under the constraint 

Having written the Eilenberger equation (|3"B")) and its 
action (j4"0)) we now use the standard technique [1^, H(| 



to derive the quantum kinetic equations. We proceed 
at one-loop order with respect to the effective interac- 
tion, which is equivalent to the "dirty RPA" [12]. The 
one-loop derivation is controlled by the small parameters 
1/Tt<£ <C 1 and a <C 1. More specifically, following the 
framework of Ref. [32| , we introduce three different distri- 
bution functions for each \i. The first one, f ll (e,x,t) 1 de- 
scribes the bare electrons, moving with velocity vf- The 
other two, Nj}(w, x, t) and N^(u), x,t), describe two types 
of bosons, having velocities u p — vf/K and u g = vf- 
The bosons of the first kind represent the usual plas- 
mons (p) of the LL, whereas those of the second kind 
are "ghosts" (g) constructed from the bare electron-hole 
pairs, thus preventing from a double-counting of the de- 
grees of freedom in the system [see the discussion around 
Eqs. @2)-(SS) below]. 

We first apply the gauge transformation 

g*{xMM) = e-^^^g»{x,t u t 2 y § " {xM) , (42) 



j* + (T x 8 2 has the same Keldysh structure 



where 9 11 

as the field <fi. This transformation is similar to that in 
Ref. [Ii,[iH, but different in that the equation of motion 
for the phase 9 in the field <j> will incorporate disorder, see 
Eq. (|5Tj) below. The "rotated" Green's functions g^ are 
expressed in terms of the electron distributions f^(x,t), 
written in the time domain, as 



5(h-t 2 ) 2W t (t 1 ,t 2 ,x) 
S(h-t 2 ) 

where h» = <% - t 2 ) - 2f(t 1 ,t 2 ,x), 
de 



2tt 



and we impose the condition 
f^(ti,h,x)\ tl ^ t2 = 



e«(*i-*»)//'[ a!> (i 1 +f 2 )/2] 



27r(ti - 1 2 + i0) 



(43) 



(44) 



(45) 



The fast charge and current fluctuations are now encoded 
in the fluctuations of the phase factors e in Eq. (|42[) . 
The gauge-transformed action reads 



s[e,j>,g] 



— S e + Sf, + Oj n t + S lmp , 



-— Tr [id t -L 6-(f> )f z g 
2vf 



Tr g T d x T 



Shit 

^iinn 



2itvf 



-Tr 



-Tr 



-(d t 0) L a x 9 + 4> a x d f 1 



f " x 2hvf 



8vpT 



-Tr e 



-i{0- 



(46) 



(47) 



(48) 
(49) 



(50) 
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where Lq = dt + f z vpd x . 

We treat the fluctuations of 9 and </> in the Gaussian 
approximation by expanding Eq. around the saddle 
point of S for a given gf. Optimizing then the action 
with respect to 9 for a given </>, we get a linear relation 
between 9 and 6: 



convenient to define the correlator of the phase $ as a 
difference of two terms: 



where 



(59) 



D~ Y 9 = -a x(i 



(51) 



where we introduce the vector notation 9 — 
(9+,9^,9+,9z) T , <j> = (</>+, 0r,^,0 2 -) T , T stands for 
transposition, and the particle-hole propagator D g is con- 
structed as 



V 



(d t + r z v F d x )a x + ^7(1 



with 



1 
r 
1 



-1 

1 2B„ 



B u — - — 



2u ^ 



de (1 - b-* u 



(52) 

(53) 
(54) 



The scattering operator 7 (|53")) describes the de- 
cay/recombination of the collective electron-hole excita- 
tion into/from the electron and hole moving in opposite 
directions, assisted by impurity scattering. Note that the 
app roximation (|5Tj) is equivalent [13] to the "dirty RPA" 
03 in Sec. IV. 

Substituting Eq. (f5Tj) back into the approximate 
quadratic action, we obtain the "dirty-RPA" propagator 
of the effective interaction 



T ) = ~v = - [& S T X Vf x - n 



where 



n 



2ttvf 



Ox Ox — O x 



(55) 



(56) 



is the polarization operator. By combining Eqs. ((55) and 
(ISTj) we get the correlator of the phases 9 (cf. Ref. 43] ) 



D q a x V a x D q 



ITTVp 



D„ - D, 



(57) 



where D p is the renormalized particle-hole propagator 
corresponding to the plasmon modes with velocity u 
given by Eq. ((6]): 



d: 



VFTzdz 



~7(1 ~ fx) • 



(58) 



As follows from Eq. ([50]) . the only phase that is coupled 
to the electron backscattering off disorder is $ = \(9~ — 
9 + ). Another observation is that the propagators of the 
fluctuations of 9 have two different types of poles: to = 
±uq and u> = ±i>f9, both smeared by disorder. It is thus 



lITVp 

^2d7 



(60) 



//// 



and b = p,g denotes the plasmon and ghost modes, which 
differ from each other in that the plasmon mode is char- 
acterized by velocity u, whereas the ghost mode by ve- 
locity vp- Then the Wigner-transform of the Keldysh 
correlator ($$ t )k has to be described by four different 
distribution function, N^(ui, x, t) and N^(uj, x, t), evolv- 
ing with velocities u and to the right and to the left: 

{tt T ) K (w,q~±u/u,X,t) = (61) 

(<S><S> T ) K (Lj,q~±iu/vp,x,t) = (62) 
- [2N±(w,x,t) + l\ TmCf(u,q) ■ 

To derive the kinetic equation for the electron distribu- 
tion function, the next step is to write down the equation 
of motion for the gauge-transformed Green's function <p 
(1421). The latter follows from the relation 



(S e + "Simp) — 

)nf J - 



Using Eq. (f5l"j) we represent S e as 



S e = _ ^7- Tr (idtTz - <f>J g 
where the shifted phase (j) = + 4>2 o x , 

4>a = ^2(Oxj)aP$l3 , 



(63) 



(64) 



(65) 



and 7 is given by Eq. (|53"|) . Notice that the field <f> does 
not depend on the chiral index /i. The Eilenberger equa- 
tion for g^ thus reads 



ip,v F d x g^ 



idt-^+^e^g-^e- 2 ^, 
4r 



.9" 



= 



(66) 

Equation (f6"6")) has to be averaged over the fluctuations of 
the phase $ with the correlator given by Eq. (|59[) . Within 
the "dirty-RPA" it is sufficient to represent g^ as a sum 
9^ = (<? M ) + where (<? M ) is the mean value, and take 
into account only the term in 5g^ that is linear in the 
fluctuations of $, keeping in mind that the quadratic- in- 
$ fluctuations of <p are incorporated in the mean value. 

By linearizing Eq. (I66p around the average (g^ 1 ), we 
obtain 8g^ — — 2(5/^) cr + , where the fluctuation of the 
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distribution function obeys 



The elastic (antisymmetric in /i) part is given by 



E 



SF(e u e 2 ) = — ^ 

0=1,2 



A£(ei,e 2 ) 

(67) 

In Eq. (|6"7| , <5/ M (ei,62) is the Fourier transform of 
5f^{t\,t2,x) and = t\ — e 2 . The source terms A^ are 
expressed through the averages as 



A?(ci,e 2 ) 
A2(ei,e 2 ) 



m(^ 2 -<) 



-e 



1 + (^-^J 

^AT-^A7 



(68) 
(69) 



where h» = 1 - 2/f . 

Notice that the general formalism of the "nonequilib- 
rium functional bosonization" [Eq. (|66p ] allows, in prin- 
ciple, for a nonperturbative treatment of both the elastic 
renormalization and the inelastic scattering if the phases 
3> are kept in the exponents (see, in particular, the calcu- 
lation of the tunneling density of states in the clean LL 
out of equilibrium in Ref. (34[). For our purposes in this 
paper, it is sufficient to expand the exponential factors 
to second order in <f>. 

The Eilenberger-type equation for the average (g 1 *) can 
now be written in the form 



id. 



4r 



= st: 



st: 



other, St e _ e , stems from the contraction of the linear 
correction dg^ with the fluctuations of <F 

The kinetic equation for the electron distribution func- 
tion is obtained by taking the Keldysh part of Eq. (|70l) : 



[d t + fiv F (d x 

+ St e _ e 4 



E 

b=p,g 



-— (7" 

+ StM 



(71) 



The electron-boson collision integral Stg_ h describes 
emission and absorption of the bosons of type b = p, g 
by the fermions. Its inelastic part, which is symmetric 
with respect to the channel index /i, reads 

SC_!b(e) = ijE / duj u)]C e -b{uj) 

x [NZMfZMl ~ f!?) - [1 + Nj!(u>)]f? (1 - /.-"J] • 

(72) 



Stl b (e) 



± 



4^ 



(73) 



where the signs ± refer to the plasmons (+) and ghosts 
(— ). In Eqs. ([7T|) - (f73|) . r and u are understood as renor- 
malized down to a scale wo, so chosen that Ae <C ujq <C A 
but a In (wo /Ae) <C 1, where Ae is a characteristic energy 
scale for the distribution function (for the quantum wire 
biased by a voltage V it is given by max{T, eV}). The 
high energy renormalization (see Sec. Ill), coming from 
scales larger than loq and independent of the details of the 
nonequilibrium state at low energies, can thus be taken 
into account at the level of input parameters (r and it) 
for the kinetic equations. 

The e-e collision integral St e _ e describes inelastic 
fcrmionic collisions due to the interaction via the plas- 
mon waves: 



-b ■ 

(70) 

The collision integrals in the right-hand side of Eq. (|7D|) 
come from the averages of second order in the fluctua- 
tions of There are two types of the collision integrals. 
One, St e _ b , comes from the second-order terms in the 

expansion of the phase factors e ±2lM * in Eq. (|66p. The where 



St e _ e (e) = tE / dude' K, e - e {uj) 

x [/ra-^-j/e-ta-zr) 
- ft-M - CT/fa - fr-H)] • (74) 

Note that St e - e does not depend on the chiral indices. 
The collision kernels are written as 

K e . b {u) = ±^f [ p-ReJ2^ v (u;, q ) , (75) 



IC e - e (uj) = K{u>) + K. e ~ g (uj) - /C e _p(w) 



(76) 



^M = ~ / ^E Re ^n^)Im^>,g) (77) 

and the explicit form of the propagators D^ u and V^ v 
can be found in Ref. [see Eqs. (A16),(A17),(A21)- 
(A23) there; note that D£" corresponds to vpD^ in 
Eqs. (A16),(A17)]. 

For the electron-boson terms we obtain in the ballistic 
limit of energy transfers u> larger than t _1 the simple 
expressions 



K, e - P {uj) 



2K 



/C e _ 9 (w) = 



(78) 



The asymptotic behavior of JC(u>) in three parametrically 
different ranges of u> in the limit of weak interaction a <C 
1 is as follows: 



/C(w) 



2a 2 /uj 2 T, t- 1 < uj < aTi , 
8a 4 r, aT x < u < Tj , 
2/w 2 t, w > Ti , 



(79) 
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FIG. 2: Frequency dependence of the collision kernel lC(u) 
[Eq. (J77J)] for q = 0.05. Three different types of scaling be- 
havior of the dimensionless product t -1 /C(u;) are indicated, 
as well as the characteristic values of u . 



where T\ ~ l/a 2 T is the characteristic temperature be- 
low which the localization effects are strong (see the dis- 
cussion at the end of Sec. IV). The log-log plot of K,(oj) 
in the whole range of u> for a particular value of a (taken 
very small for the purpose of illustration) is shown in 

ng.rj 

An important feature of K,(ui) in Eq. l|79[) is its nonper- 
turbative behavior with respect to a at w > aT\ ~ 1/cvt. 
In particular, for uj » T\ the e-e collision kernel does 
not at all depend on the e-e interaction strength. The 
origin of the nonperturbative dependence on a, despite 
q« 1 being small, is related to the analytical structure 
of K(oj) in Eq. (|77|) . Specifically, the integrand of JC(u) 
contains eight poles, which in the limit of large w are only 
slightly "damped" by disorder: q ~ ±w(l ± i/2u>r)/vF 
and q ~ ±w(l ± i/2ujt)/u. As a result, the contour of in- 
tegration in the plane of q is squeezed between two close 
poles (u — * vp for a — ► 0), one of which is in the upper- 
half plane and the other in the lower one. At uj ^> T\ we 
thus have K.(u>) oc a 2 /\u — Vf\ and a drops out altogether. 

The kinetic equations for the bosonic distribution func- 
tions N£ g follow from Eqs. ([59]) . (I6T|) . (j62|) : 

(d t + fiu b d x )Ng(u) = --N£(w) + St 6 _ e (u;) , (80) 

T 

where 

St6-eH = ^ £ / ~ /e-t) (81) 

describes the creation of the boson from an electron-hole 
pair, where the electron and hole move in the opposite 
directions. 

The role of the ghost modes can be further elucidated 
if one considers the energy conservation law. The elec- 
tronic and bosonic energy densities p\, p\ and the energy 



current densities j|, j% are given by 

l r°° 

2n 



Je 

p\ 



1 f°° 

s J 



1 



3b 



27T-U& J 

1 

2^ 



dojw[N+(oj) + N b -(uj)] 
da;w[JV^(«)-iV 6 -( W )] . 



(82) 
(83) 
(84) 
(85) 



The prefactors in Eqs. (|82j) and (|84|) represent the density 
of states for electrons and bosons, respectively. Then the 
kinetic equations ([7T|) and (f80|) assure the conservation 
law 



(86) 



" Pp ~ Pg) + 9 x (fe + fp ~ fg) =je£, 

where in the right-hand side £ is the applied electric field 
and j e is the induced current of charge. The kinetic en- 
ergy without any interaction is determined by the energy 
of the bare electrons. In the presence of Coulomb in- 
teraction, the plasmon energy is given by a sum of the 
e-e interaction energy and the kinetic energy of the bare 
electron-hole pairs, the latter being the ghost energy by 
construction. The total energy is thus given by the sum 
of the plasmon and electron systems with a subtraction 
of the energy of the ghosts. 

We now turn to the rate of energy relaxation t^ 1 in 
the limit of weak nonequilibrium by linearizing the ki- 
netic equations (|7ip and (|5Tj) . One sees that at large en- 
ergy transfers the inelastic e-e scattering dominates over 
electron-boson collisions: K(uj) S> JC e - g {u) — JC e ~b(uj) 
for w l/a z / 2 T. Assuming that the large uj give 
the main contribution to the energy relaxation in the 
limit T > Ti, where the localization effects can be ne- 
glected (see Sec. IV), the equilibration rate at which the 
fcrmionic system relaxes to a locally equilibrium Fermi 
distribution is estimated as 

rT , 

dujuj 2 K{uj) ~ T 2 K(T) ~ - . (87) 



1 



t e {t) : 

Notice that the characteristic uj in Eq. (|87|) is of order 
T, which justifies the use of fC(u>) only. On the other 
hand, this makes it impossible to describe the equilibra- 
tion in terms of the much simpler Fokker-Planck equation 
in the energy space. Remarkably, the equilibration rate 
(|87p does not depend on the strength of interaction and 
is given by the elastic scattering rate. The interaction 
constant a enters only through the condition T 3> T±. 
The equilibration rate turns out to be much smaller than 
the (clean) e-e collision rate (f22|) and also much smaller 
than the dephasing rate (|3T))) : 

r E x «r - 1 «r- 1 . (88) 

It is also worth emphasizing that the characteristic en- 
ergy transfers are parametrically different in these three 
types of relaxation processes. 
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VI. SUMMARY 

In this paper, we have studied the relaxation proper- 
ties of interacting spinless electrons in a disordered quan- 
tum wire within the Luttinger-liquid model. We first re- 
view the basic concepts in the disordered Luttinger liquid 
at equilibrium, including the elastic renormalization, de- 
phasing, and interference-induced localization. We have 
introduced the general framework for describing the re- 
laxation processes in the strongly correlated (non-Fermi) 
electron system at nonequilibrium. Our main result is the 
coupled set of the kinetic equations for the fermionic (|TTj) 
and bosonic ([50)1 distribution functions. The peculiarity 
of the Luttinger liquid model is that the electron-electron 
scattering rate (the inverse lifetime of the fermionic ex- 
citations) is finite at nonzero temperature but the en- 
ergy exchange is exactly zero in the clean limit. The 
energy relaxation occurs only due to the scattering off 
disorder. We have calculated the energy equilibration 
rate that turns out to be independent of the strength 
of electron-electron interaction at sufficiently high tem- 
perature, when the Anderson localization effects are sup- 
pressed, and equal to the rate of elastic scattering off 
disorder. 
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